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This  report  reviews  the  existing  state  of  knowledge  concern¬ 
ing  the  evaluation  of  the  forces  imposed  on  a  body  in  a  fluid 
due  to  acceleration  of  either  the  body  or  the  fluid.  It 
concentrates  on  those  fluid  inertial  forces  due  to  acceleration 
rather  than  on  the  drag/lift  forces  due  to  steady  motion.  The 
first  part  of  the  report  presents  a  survey  of  the  analytical 
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geometry  are  explored  in  a  preliminary  way.,  The  possible 
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further  highlighted.  The  confused  state  of  affairs  regarding 
the  possibly  major  effects  of  viscosity  in  certain  regimes  of 
frequency  and  Reynolds  number  is  discussed.  Finally  a  number 
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put  forward. 
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L  INTRODUCTION 

W  henever  acceleration  is  imposed  on  a  Quid  flow  either  by  acceleration  of  a  body  or  by  acceleration 
externally  imposed  on  the  fluid  additional  fluid  forces  will  act  on  the  surfaces  in  mntart  with  the  fluid 
These  fluid  inertial  forces  can  be  of  considerable  importance  in  many  ocean  engineering  problems.  The 
purpose  of  this  report  is  to  review  some  of  the  characteristics  of  these  fluid  inertial  forces  and  in  partic¬ 
ular,  to  evaluate  the  state  of  knowledge  of  the  "added  mass”  mainoes  which  are  used  to  characterize  the 
forces  The  first  part  of  the  report  (Section  3)  is  also  intended  to  serve  educational  purposes.  The 
second  part  (Section  4)  reviews  the  existing  data  base  end  some  of  the  areas  in  which  there  is  either  a 
lack  of  data  or  a  data  base  which  is  contradictory.  It  is  also  intended  to  convey  the  limitjetnnn  of  the 
existing  knowledge.  Finally  a  number  of  suggestions  for  improvement  in  our  present  understanding  are 
listed  in  the  concluding  section 

Unlike  many  reviews,  the  author  has  not  attempted  to  absorb  every  publication  on  the  subject. 
Rather  the  time  which  would  have  been  spent  on  such  an  effort,  was  devoted  to  more  concentrated 
analysts  of  the  subject.  Other  excellent  reviews  of  vanous  aspects  of  unsteady  fluid  forces  exist;  in  par¬ 
ticular  the  reader  is  referred  to  the  recent  books  by  Blevins  (Ref. IB)  and  by  Sarpkaya  and  Isaacson 
(Ref.  17). 


Z  GENERAL  EXPLANATION  or  ADDED  BASS 


Perhaps  the  simplest  view  of  the  phenomenon  of  "added  mans"  is  that  it  determines  the  necessary  work 
done  to  change  the  kinetic  energy  associated  with  the  motion  of  the  fluid.  Any  motion  of  a  fluid  such 
as  that  which  occurs  when  a  body  moves  through  it  implies  a  certain  positive,  non- zero  amount  of 
kinetic  energy  associated  with  the  fluid  motions.  This  kinetic  energy,  T,  can  be  amply  represented  by 

T=  +u£)dV=  | -jfmndV  ( 1) 

where  the  Ui(i=l,  2, 3)  represent  the  Cartesian  components  of  fluid  velocity  and  V  is  the  entire 
domain  or  volume  of  fluid.  For  simplicity  we  shall  assume  throughout  that  the  fluid  is  incompressible 
with  a  denaty  p. 

If  the  motion  of  the  body  is  one  of  steady  rectilinear  translation  at  velocity.  U.  through  a  fluid  oth¬ 
erwise  at  rest  then  dearly  the  amount  of  kinetic  energy,  T,  remains  constant  with  time.  Furthermore  it 
is  deer  that  T  will  in  some  manner  be  proportional  to  the  square  of  the  velodty,  U,  of  translation 
Indeed  if  the  flow  is  such  that  when  U  is  altered  the  velocity,  14.  at  each  point  in  fluid  relative  to  the 
body  varies  in  direct  proportion  to  U  then  T  could  conveniently  be  expressed  as 

T-pjtfl  uhem  (2) 

and  the  integral  /  would  be  a  simple  invariant  number,  This  is  indeed  the  case  with  some  fluid  flow 
solutions  such  as  potential  flow  and  low  Reynolds  number  Stokes  flow.  However  it  may  not  be  true  for 
the  complex,  vortex  shedding  flows  which  occur  at  intermediate  Reynolds  numbers. 

N  ow  consider  that  the  body  begins  to  accelerate  or  decelerate.  Clearly  the  kinetic  energy  in  the  fluid 
will  also  begin  to  change  as  U  changes.  If  the  body  is  accelerated  then  the  kinetic  energy  will  in  all 
probability  increase  But  this  energy  must  be  supplied;  additional  work  must  be  done  on  the  fluid  by 
the  body  In  order  to  increase  the  kinetic  energy  of  the  fluid.  A  nd  the  rate  of  additional  work  required 
is  simply  the  rate  of  change  of  T  with  respect  to  time,  dT  ML  This  additional  work  is  therefore  experi¬ 
enced  by  the  body  as  an  additional  drag,  F,  such  that  the  rate  of  additional  work  done,  -FV  is  simply 
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equal  to  dTAB.  It  the  pattern  of  flow  is  not  changing  such  that  the  integral  I  remains  constant  it  fol¬ 
lows  that  the  "added  dreg”,/!  is  simply 


U  dt  p  dt 


(3) 


Now  this  tone  has  the  same  form  and  agn  as  that  required  to  accelerate  the  mass  (m)  of  the  body 

itself,  namely  m^-.  Consequently  it  is  often  convenient  to  visualize  the  mass  of  fluid  pi,  as  an  "added 
at 

mass",  M,  of  fluid  which  is  being  accelerated  with  the  body.  Of  course,  there  is  no  such  identifiable 
fluid  mass;  rather  all  of  the  fluid  is  accelerating  to  some  degree  such  that  the  total  kinetic  energy  of  the 
fluid  is  increasing 


It  is  important  to  stress  that  Fis  not  the  only  drag  force  experienced  by  the  body.  During  steady 
translation  in  a  real  viscous  fluid  there  is  a  steady  dreg  associated  with  the  necessary  work  which  must 
be  done  to  balance  the  steady  rate  of  dissipation  of  energy  in  the  viscous  fluid  When  the  body 
accelerates  there  will  be  a  similar  though  not  necessarily  equal  drag  associated  with  the  instantaneous 
value  of  U.  Furthermore  there  may  be  delayed  "fleets  associated  with  the  entire  previous  history  of 
trandation  (e  g.  the  Basset  force,  Ref.l,  p.375) 


-  ■.  aUiMI 
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a  ANALYTICAL  APPROACHES  TO  ADDED  HASS 

at  EXAMPLES:  RECTILINEAR  MOTION  OPA  SPH HIE  AND  CYLINDER  WITH  POTENTIAL  PLOW 


In  the  preceding  discussion  the  consequences  of  acceleration  were  illustrated  by  reference  to  simple  rec¬ 
tilinear  motion  of  velocity,  U.  It  should  be  dear  that  the  methodology  could  be  extended  to  more  gen¬ 
eral  motions  and  indeed  this  will  be  carried  out  in  the  following  section.  But  prior  to  this  it  is  worth 
illustrating  how  the  integral,  /,  and  therefore  the  added  mass  can  be  calculated  for  rectilinear  motion 
For  the  purposes  of  this  example  let  us  examine  the  idealized  potential  flows  pest  a  sphere  and  a 
cylinder  The  geometry  for  both  is  as  depicted  in  Fig.l.  The  sphere  or  cylinder  of  radius  R  is  assumed 
to  be  moving  with  time  varying  velocity  U(t)  (I  is  time)  in  the  positive  x  direction.  Polar  coordi¬ 
nates  (r,i$)  are  used  where  z=rcosiS 


Fig. 


The  resulting  fluid  velocities  Ur/u*  in  the  r  and  i?  directions  are  then  given  by  a  velocity  potential?  , 
such  that 


3?  1  dtp 

^=77  :  ^"7  ™ 


(4) 


and  the  appropriate  velocity  potentials  in  the  two  cases  are 

...  UR3 


2r* 


oos 


(5) 


?t^dr=--^-cosi> 


The  reader  who  is  unfamiliar  with  these  solutions  may  wish  to  satisfy  himself  that  two  solutions  satisfy 
(i)  Laplaces  equation  V  8?=0,  in  spherical  and  cylindrical  coordinates  respectively  and  (ii)  the  boun¬ 
dary  condition  that  the  relative  velocity  normal  to  the  surface  of  the  body  is  zero  ((14 )rmtf=U  oos 


-5- 


It  follows  that  these  flows  are  of  the  type  in  which  14  is  directly  proportional  to  U  and  consequently 
the  integrals,  /,  can  be  evaluated  as 

^J=ff  FIT  +  iw  Srr^sin d&dr- rr/?3  (7) 

therefore  M  =  rr/?3p 

" Zn  1  1Z  f  1  fl  1*1 

Ctprder  I=f  f  ~rr  Trfi*fr=rr.ff2  per  lent  length  (B) 

therefore  M  =nRzp  per  writ  Length. 

Note  that  the  added  mass,  M,  of  the  cylinder  is  equal  to  the  mass  of  the  fluid  displaced  by  the  body, 
whereas  the  added  mass  of  the  sphere  is  one  half  of  the  displaced  mass. 

These  are  algebraically  the  simplest  potential  flows  for  which  the  value  of  /  can  be  evaluated  How¬ 
ever.  it  is  conceptually  ample  to  visualize  how  /  ^uld  be  evaluated  for  any  flow  provided  the  flow  solu¬ 
tion  (the  Ui/U  values)  are  available.  Note  that,  in  effect  one  need  only  have  available  the  solution  for 
the  steady  flow  in  the  direction  under  consideration.  This  considerably  amplifies  the  added  mass  calcu¬ 
lation  for  rectilinear  motion.  Later  we  shall  examine  the  more  general  case  of  arbitrary  motion 

a2  RELATION  TO  DISPLACED  MASS;  VARIATION  WITH  DIRECTION  OF  ACCELERATION 

In  the  preceding  section  it  was  noted  that  in  the  ideal  case  of  potential  flow  around  a  droular  cylinder  in 
rectilinear  motion  the  added  mass  is  equal  to  the  mass  of  fluid  displaced  by  the  cylinder.  This  should 
be  regarded  as  coincidental.  There  is  no  general  correlation  between  added  mass  and  displaced  fluid 
mass.  As  we  have  seen  the  added  mass  for  a  sphere  is  one  half  of  the  displaced  fluid  mass.  Further¬ 
more  the  idealized  potential  flow  past  an  infinitely  thin  flat  plate  (zero  displaced  fluid  mass)  accelerated 
normal  to  itself  has  an  added  mass  equal  to  the  mass  of  a  circular  cylinder  of  fluid  with  a  diameter 
equal  to  the  width  of  the  plate. 

Thus  the  displaced  fluid  mass  may  not  even  be  a  good  first  approximation  to  the  added  mass  (except 
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for  the  case  of  the  circular  cylinder) .  Furthermore  we  shell  see  that,  in  general,  the  value  of  the  added 
mass  depends  on  the  direction  of  acceleration  For  example,  the  idealized  potential  flow  solution  for 
the  infinitely  thin  flat  plate  accelerated  in  a  tangential  rather  than  a  normal  direction  yields  zero  added 
mass  rather  than  the  value  described  above.  A  review  of  the  available  date  suggests  that  a  better  (but 
still  very  crude)  first  approximation  to  the  added  mass  of  a  body  for  a  given  direction  of  acceleration 
would  be  the  mass  of  the  fluid  volume  obtained  by  taking  the  projected  area  of  the  body  in  that  direc¬ 
tion  and  evaluating  one  half  of  the  volume  of  the  sphere  with  the  same  projected  area  (see  Sections 
4.2,  4.3).  An  improvement  on  this  is  included  in  Section  4.2. 

One  other  complication  will  emerge  in  the  following  section  when  the  complete  added  mass  matrix 
is  defined,  namely  that  the  force  on  the  body  due  to  acceleration  is  not  necessarily  in  the  same  direction 
as  the  acceleration.  For  an  unsymmetric  body  acceleration  in  one  direction  can  give  rise  to  an  "added 
mass'  effect  resulting  in  a  force  which  has  a  component  in  a  direction  perpendicular  to  the  direction  of 
acceleration.  If,  for  example,  one  were  lifting  a  body  from  the  ocean  bottom  by  means  of  a  cable  then 
an  increase  in  the  lift  rate  could  produce  a  lateral  otion  of  the  body. 

a3  THE  ADDED  HASS  MATRIX 

Up  to  this  point,  most  of  the  examples  and  discussion  have  centered  on  simple  rectlinear  motion 
However  in  general  the  response  of  a  body  to  an  additional  force  applied  at  some  point  and  in  some 
direction  will  not  be  confined  to  motion  in  that  same  direction  Instead  there  will  be  a  general  induced 
acceleration  of  the  body  consisting  of  three  translation  accelerations,  Aj,  j= 1, 2, 3  in  three  perpendicu¬ 
lar  directions  and  three  angular  accelerations, A},  j- 4,5,6.  Then  the  added  mass  matrix 
A/*/,  i=l-*  6,  provides  a  method  of  expressing  the  relationship  between  the  six  force  com¬ 

ponents.  Ft.  imposed  on  the  body  by  the  inertial  effects  of  the  fluid  due  to  the  six  possible  components 
of  acceleration, : 

Fi=~M1iAi  (9) 

The  matrix  A/„  must  have  added  to  it  the  inertial  matrix  due  to  the  mass  of  the  body  in  order  to  com¬ 
plete  the  formulation  of  the  inertial  forces  If  the  center  of  mass  of  the  body  is  chosen  as  the  origin  the 
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body  mass  matrix  is  symmetric  and  contains  only  seven  different,  non-zero  values,  namely  the  masB 
and  the  six  different  components  of  the  moment-of-mertia  matrix  [Yih,  p.102].  However  one  cannot 
in  general  relate  any  of  the  36  different  components  of  the  added  ma®  matrix  nor  prove  that  any  of 
them  are  zero  except  in  specific  cases  or  for  specific  kinds  of  flow.  Consequently  an  externally  applied 
additional  force  will  in  general  create  acceleration  in  all  six  components  of  velocity  and  angular  velocity. 
Thankfully  it  is  rarely  necessary  to  have  to  hen  die  36  different  added  mass  coefficients.  For  potential 
flow  one  can  show  [Yih.  p.100]  that,  the  added  mass  matrix  must  be  symmetnc;  since  the  system  is 
then  conservative  the  symmetry  also  follows  from  the  theorem  of  reciprocity.  This  reduces  the  number 
of  coefficients  to  21.  However  no  further  reduction  is  possible  except  for  bodies  with  geometric  sym¬ 
metries 

The  simplifications  introduced  by  geometric  symmetries  of  the  body  are  fairly  easily  established. 
Consider  for  example  a  body  with  a  single  plane  of  symmetry,  for  example  an  airplane  It  is  dearly 
convenient  to  select  axes  such  that  this  plane  of  symmetry  corresponds  say,  the  x3=0  plane.  Then  any 
acceleration  confined  to  this  plane,  namely  any  -  'mbination  of  A  LAS  and  Ae  will  produce  no  added 
mass  force  F3.  FA  or  Fb\  the  only  possible  non-zero  forces  will  be  Fv  Fz  and  Fe  It  follows  that  for 
such  a  body  the  following  9  components  of  the  added  mass  matrix  will  be  zero: 

Mq=0  for  i=3,4,  5,  j- 1,2,6  (10) 

If  in  addition  the  flow  is  assumed  to  be  potential  such  that  the  matrix  is  symmetric  then  A/j,=0  for  the 
same  domains  of  i  and  J.  The  number  of  non-zero  values  required  to  define  the  matrix  is  IS  namely 

M±<  i=l-*6  and  ^12,^34,^35,^40.^16  and  M gg  (11) 

Bodies  which  have  two  planes  of  symmetry  (for  example  a  hemisphere)  yield  a  further  reduction  In 
the  number  of  non- zero  values.  Suppose  axes  are  chosen  such  that  both  xz=0  and  x3=0  ere  planes  of 
symmetry  Then  not  only  must  (10)  be  true  but  also 

Mv=0  for  i=2, 4,6  ;  ;  =  1,3, 5  (12) 


and  eg  am.  assuming  potential  flow  0  for  the  same  domains  Then  the  only  non-zero  values  which 
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need  evaluation  ere 

tf4.  i=l-8  end  M » *  (13) 

The  lest  two,  which  with  M  &=M?£  and  represent  the  only  non-zero  off-diagonal  terms, 

correspond  to  the  moment  about  the  x3  axis  generated  by  acceleration  in  the  direction  end  the 
moment  about  the  x2  axis  generated  by  acceleration  in  the  x3  direction  In  other  words  since  the 
body  is  not  symmetric  about  the  xa x3  plane  linear  acceleration  in  either  the  xa  or  x3  direction  will 
cause  pitching  moments  in  the  x,x2  or  x,x3  planes. 

A  few  simple  bodies  such  as  a  sphere,  circular  cylinder,  cube,  rectangular  box,  etc  hare  three  planes 
of  symmetry  By  following  the  same  procedure  used  above  it  is  dear  that  the  only  possible  non-zero 
elements  are 

My,  i=l-*6,  U  is, A#  ifl, -W  gf.A/gg.Af  31,/Zgg  (  14) 

and  that  if  potential  flow  is  assumed  all  of  the  off-diagonal  terms  are  zero.  Only  in  this  ample  case  of 
three  axes  of  symmetry  and  symmetry  of  the  nrv'.  x  (see  below  for  conditions  on  this)  does  the  matrix 
become  purely  diagonal  so  that  there  are  no  secondary  induced  accelerations. 

Jt  remains  to  discuss  the  precise  flow  conditions  under  which  tl«.  matrix  can  be  assumed  to  be  sym¬ 
metric  and  then  finally  to  indicate  how  all  of  the  elements  could  be  evaluated 

ADDED  MASS  MATRIX  SYMMETRY  AND  SUPERPOSIBILITY  0 Y  PLOW  SOLUTIONS 

The  astute  reader  will  have  recognized  that  the  mere  definition  of  My  in  Eq.(0)  requires  certain 
assumptions  concerning  the  nature  of  the  flow  and  the  ability  to  linearly  superpose  the  effects  (i.e. 
forces)  of  acceleration  m  the  six  directions.  The  question  of  the  minimum  preconditions  necessary  in 
order  to  write  Equation  (9)  is  one  which  will  not  be  addressed  here.  It  is  however  deer  that,  these 
preconditions  are  met  as  soon  as  one  makes  the  assumptions  necessary  to  evaluate  My  To  the  authors 
knowledge  the  only  evaluations  which  exist  require  that  the  fluid  flow  is  superposable  in  the  following 
sense;  that  the  total  induced  fluid  velocity  can  be  obtained  by  linear  addition  of  the  fluid  velocities 
caused  by  each  of  the  components  of  the  body  motion  or  velocity.  For  this  to  be  true  requires  that 


both  the  equations  used  to  solve  for  the  fluid  flow  and  the  boundary  conditions  be  linear.  This  is  not 
true  in  general  of  the  Navier- Stokes  equations  for  fluid  motion  and  therefore  superposability  is  not  in 
general,  applicable.  However  there  are  two  models  of  fluid  flow  which  do  satisfy  this  condition  namely 
(i)  the  potential  flow  model  for  high  Reynolds  flow  [Yih,  p.100]  and  (ii)  the  Stokes  flow  model  for 
asymptotically  small  Reynolds  numbers.  In  both  cases  the  equations  of  motion  can  be  put  in  linear 
form  Furthermore  provided  if  one  is  dealing  with  rigid  or  undeformable  boundaries  the  boundary  con¬ 
ditions  are  also  linear  Only  in  these  two  limiting  oases  can  the  added  mass  matrix  be  regarded  as  an 
exact  representation  of  the  relation  between  fluid  inertial  force  and  body  acceleration.  In  other  types  of 
flow  it  oould  however  be  regarded  as  a  reasonable  first  approximation.  Case  (ii)  above  is  of  interest  in 
flows  such  as  occur  in  slurries  or  suspensions;  however  we  shall  from  here  on  confine  our  remarks  to 
case  (i)  which  is  of  greater  practical  importance  in  ocean  engineering. 

W  hen  the  flow  is  linearly  superposable,  it  is  convenient  to  define  Uy  as  the  induced  fluid  velocity 
caused  by  unit  velocity  of  the  body  in  the  j  direction  (>=  1  -*  0).  lnduded  here  are  both  the  transla¬ 
tion  components,  j  =  1 . 2. 3.  and  rotational  components,  >=4,5,6  of  body  motion  Then  if  the  body 
velocities  are  denoted  by  Uy  j- 1  -» 6.  it  follows  that  the  fluid  velocity 

Ui=y^Uj  (15) 

Consequently  one  can  wnte  Equation  (1)  as 

T=jA*UiUk  (16) 

where  the  matrix  is  composed  of  elements 

Af'=pfvqUledV=Mfe  (17) 

It  con  be  shown  [Yih,  p.  102]  that  the  matrix  A ^  is  in  fact  the  added  mass  matrix  Mp.  It  is  certainly 
dear  that  the  diagonal  terms  Axl  A&  Ak  are  identical  to  the  added  masses  evaluated  in  Section  3.  (To 
establish  this  define  the  direction  x  of  Section  3  as  either  i,  i2  or  za;  then  and  ifc  are  identical 

IL 

and  equal  to  the  veiodty  —  used  in  Section  3.) 


Furthermore  it  la  deer  from  this  evaluation  that  the  added  mass  matrix  oust  be  symmetric  Maoe 
reversing  j  and  k  in  Equation  ( 16)  does  not  change  the  value  of  the  integral.  Hence  superpoanbility 
implies  symmetry  of  the  added  mass  matrix. 

a5  EVALUATION  OF  THE  ADDED  HASS  MATRIX 

The  expression  (17)  will  permit  the  evaluation  of  the  entire  added  mass  matrix  Indeed  It  stould  be 
particularly  noted  that  use  of  this  result  only  requires  the  solution  of  steady  flow  problems  since  Uq  Is 
the  fluid  velocity  due  to  unit  velocity  of  motion  of  the  body  in  the  j  direction.  Consequently  the  solu¬ 
tion  of  six  steady  flows  for  j  =  l-»  6  allows  evaluation  of  all  21  distinct  values  in  the  added  maSB  matrix 
Hence  one  can  make  use  of  existing  methods  for  solving  steady  flows  around  bodies  of  quite  complex 
geometry  in  order  to  evaluate  the  added  mass  matrix.  References  1.2, 3, 4  and  9  provide  information  on 
these  existing  methods 

One  other  form  of  Equation  (17)  can  also  be  valuable  in  dealing  with  potential  flows.  Then  if 
represents  the  velocity  potential  of  the  steady  flow  due  to  unit  motion  of  the  body  in  the  j-diredion 
then  it  follows  that 


Substitution  into  Equation  (17)  and  application  of  Green's  theorem  leads  to 


(IB) 


(19) 


where  n  is  the  outward  normal  to  the  surface,  S,  which  represents  the  body  surface.  In  many  of 
steady  potential  flows  around  complex  bodies  it  is  dearly  easier  to  evaluate  the  surface  integral  in  ( 19) 
than  the  volume  integral  in  (17).  Indeed  the  form  (19)  is  ideally  suited  for  use  with  potential  flow 
codes  such  as  the  D  ouglas-N  euman  code. 

aa  VELOCITY  AND  ACCELERATION  OF  THE  FLUID  RATHER  THAN  THE  BODY 


A 1!  of  the  preceding  discusaon  was  centered  on  the  inertial  fortes  due  to  acceleration  of  a  body  in  a 
fluid  This  review  would  be  incomplete  without  some  comment  on  the  coses  in  which  thB  HuM  far 
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from  the  body  is  either  (i)  moving  with  a  constant,  uniform  velocity  or  (ii)  accelerating. 

Examine  case  (i)  first  It  was  implicitly  assumed  in  all  the  preceding  sections  that,  the  fluid  far  from 
the  body  was  at  rest  Otherwise  dearly  the  integral  defining  T  (Eq.(  1))  would  have  an  infinite  value 
and  the  subsequent  analysis  would  be  meaningless  If,  as  in  case  (i),  the  fluid  far  from  the  body  has 
some  uniform  constant  velocity  denoted  by  Wt  then  it  is  dear  that  since  the  inertial  force  cannot  be 
altered  by  a  simple  Galilean  transformation  it  follows  that  the  proper  definition  of  T  under  these  cir¬ 
cumstances  is 

T=%-f(ul-Wi)(ul-Wi)dV  (20) 

a  v 

The  value  of  this  integral  is  then  finite  and  the  conundrum  resolved.  In  other  words  the  appropriate  14 
to  be  used  in  Eq.(l)  is  the  velocity  of  the  fluid  •niatim  to  the  fluid  velocity  Jar  from  the  boty  provided  the 
latter  is  constant  with  time.  This  leads  to  no  alteration  in  fluid  inertial  forces.  A  rigorous  expression 
for  the  forces  would  be 

(21) 

but  since  the  time  denvative  of  Vf }  is  zero  the  original  relation  (9)  »•=  recovered. 

However  case  (ii)  in  which  W}  is  a  function  of  time  is  more  complex.  It  is  important  to  identify 
the  fluid  inertial  forces  in  this  case  for  two  reasons  First  it  is  of  practical  importance  in  analyzing,  for 
example,  ocean  wa we  forces  on  structures.  Secondly,  many  of  the  important  experiments  on  unsteady 
forces  ere  performed  uang  an  accelerating  fluid  rather  than  an  accelerating  body  (e.g.  Refs.10  and  11). 
W  e  begin  by  visualizing  a  case  (1)  flow  with  a  constant  uniform  fluid  velocity,  W j  (j=l-»  3).  far 
from  a  body  whose  center  of  volume  is  moving  at  a  velocity,  Uj  0=1-*  8).  The  body  is  also 
accelerating  with  components,  A}  The  flow  satisfies  the  Navien  Stokes  equations  for  fluid  motion  * 
and  the  solid  bocfy  boundary  conditions.  The  fluid  inertial  forces  in  this  case  are  given  by  Eq.(2l). 
Now  consider  a  slightly  different  flow  whose  velocities  are  identical  to  those  of  the  first  flow  but  in 
which  an  additional  uniform  acceleration  )  =  1, 2, 3  is  applied  globally  to  both  the  fluid  end  body.  N  ow 


the  actual  acceleration  of  the  body  ts  {Aj  +  The  NavienStokes  equations  of  fluid  motion  and 
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the  solid  body  boundary  conditions  are  identical  fa-  the  tiro  flam  im aft  that  where  the  proasure,  p 

dW, 

appears  in  the  equations  for  the  first  flow,  the  expression  p-px*  appears  in  the  equations  for  the 

second  flow.  Consequently  the  stresses  and  forces  which  the  fluid  exerts  on  the  body  ere  identical 

dWt 

except  for  an  additional  contribution  in  the  second  flow  due  to  the  additional  pressure,  px,  — — W  hen 
this  is  integrated  over  the  surface  of  the  body  the  additional  force  on  the  body  tums  out  to  be 

dWt 

p  VD  where  VD  is  the  volume  of  fluid  displaced  by  the  body.  Consequently  the  inertial  force  is 


dW 

Fi=-MqAj  +  pVD  -~- 


(22) 


But  as  stated  previously  the  acceleration  of  the  body  in  the  second  flow  is  now  A}  +  dWj/dt 
and  hence  in  the  case  of  the  second  flow 


dU, 


L- 


dt 


=  A<  + 


AW. 


(23) 


where  dW jAJt  is  the  acceleration  of  the  fluid  fa*  *rom  the  body  Substitution  for  A}  in  Eq.(22)  pro¬ 
duces  the  final  required  result  for  the  second  flow: 

dUi  ,  dWt 

Fi  =  -My  +  (My  +  p  VDSy)  — .  j=  1. 2. 3  (24) 

where  6y  is  the  Kronecker  delta  (<5^-1  for  i=j,  6y~0  for  i*j). 

Therefore  the  "added  mass"  associated  with  the  food  acceleration,  dW j/tit,  in  the  second  flow  is  the 
sum  of  the  true  added  mass  My.  and  a  diagonal  matrix  with  components  equal  to  the  mass  of  the  dis¬ 
placed  fluid,  pVD. 

However  we  must  now  examine  more  dosely  the  general  validity  of  Eq.(34).  The  first  and  second 
flows  described  above  were  assumed  to  have  identical  fluid  velocity  fields  at  the  moment  at  which  the 
forces  were  considered  This  will  ret  be  true  in  general  for  solutions  of  the  Navier-Stokes  equation 
even  though  the  body  velocities  and  far  field  fluid  velocities  are  identical.  In  general  the  solutions  to 
the  Navier-Stokes  equations  will  also  depend  on  all  of  the  previous  time  history  of  the  body  arxi  far 
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field  fluid  motions  and  consequently  the  two  flows  will  ref  in  general  have  identical  fluid  velocity  fields 
There  are  however  two  important  exceptions  to  this  and  in  both  oases  Eq.(24)  will  be  true.  First  if  the 
viscous  effects  are  neglected  then  the  fluid  has  no  memory  and  the  fluid  velocity  Adds  will  be  identical; 
thus  Eq.(24)  holds  for  potential  flows.  Secondly  if  the  previous  history  of  the  ratafee  velocity, 

( Uj  -  Wj)  is  identical  in  the  two  flows  then  (24)  will  hold  regardless  of  viscous  effects. 

Therefore,  in  summary,  the  fluid  inertial  forces  due  to  any  combination  of  body  or  far  field  fluid 
acceleration  ( dUj/dt  or  dW }/dt)  can  be  exactly  represented  by  Eq.(24)  if  either  (i)  viscous  effects 
am  neglected  or  (li)  the  matrix  Af,  represents  the  fluid  inertial  forces  for  the  case  in  which  the  fluid  is 
at  rest  far  from  the  body  and  the  entire  previous  history  of  the  relative  motion  ( U}  -  is  identical 
to  that  of  the  flow  under  oonsderation  The  latter  is  indeed  the  case  when  comparing  two  cases,  for 
example,  in  the  first  of  which  the  far  field  fluid  motion  is  sinusoidal  in  time  and  the  body  at  rest  and  in 
the  second  of  which  the  far  field  fluid  motion  is  at  rest  and  the  body  moves  sinusoidally  Consequently 
the  "added  mass'"  in  the  experiments  of  Keulegan  and  Carpenter  (Ref.  10)  in  which  the  far  field  fluid  is 
accelerated  sinusoidally  should  yield  (Mu  +  pH)  whereas  the  experiments  of  Stop,  Romberg  and 
Ferer  (Ref.  15)  in  which  the  body  is  accelerated  should  yield.  Mu  To  transfer  results  from  one  case  to 
the  other  requires  the  addition  or  subtraction  of  the  displaced  ma~:'  For  the  examples  of  Section  3.1 
the  values  of  (Af*  +  p  Vd)  would  be  2pnRz  per  unit  length  in  the  case  of  the  cylinder  and  ZpnR* 
in  the  case  of  the  sphere.  Sometimes  the  total  (Mu  +  P^d)  ls  referred  to  as  the  added  mass  and  this 
can  result  in  some  oonfusion  Strictly  speaking  the  term  added  mass  ^iould  be  reserved  for  Mq  only, 
or  in  other  words  the  case  in  which  the  body  is  accelerating  and  not  the  far  field  fluid 

a7  THE  EFFECT  OF  A  NEARLY  SOLID  BOUNDARY 

The  effects  on  the  added  mass  due  to  the  proximity  of  a  solid  boundary  will  be  addressed  in  more  detail 
later  (see  Section  4.4).  It  is  generally  true  that  the  presence  of  the  boundary  tends  to  increase  the 
added  mass  (see  Tables  I  -*  V)  and  sometimes  this  increase  can  be  very  large.  Here  we  merely  remark 
that  the  preceding  theoretical  results  are  equally  applicable  in  the  presence  of  a  solid  boundary  with  the 
following  addenda 


A.  The  reductions  due  to  geometric  symmetries  discussed  in  Section  3.3  only  apply  to  total  geometric 


synxnetries  of  both  the  body  and  solid  boundary 

B  Potential  flows  with  a  plane  solid  boundary  can  be  modelled  by  reflecting  the  flow  and  body  in  the 
plane  aid  treating  the  total  flow  due  to  the  body  and  its  im^e  Equivalence  of  the  two  problems 
allows  the  transference  of  added  mass  coefficients  from  one  to  the  other.  As  an  example  of  this 
see  the  cases  of  two  cylinders  and  a  cylinder  plus  a  plane  boundary  in  Table  II. 

a8  THE  EFFECT  OF  A  NEARLY  FREE  SURFACE 

Unlike  the  presence  of  a  solid  boundary,  a  free  surface  boundary  adds  considerably  to  the  complexity  of 
the  problem  This  is  due  to  the  fact  that  in  general,  the  boundary  condition  is  non-linear  end  hence 
superposability  is  not  satisfied  Asa  consequence  the  dynamic  behavior  of  bodies  near  a  free  surface  is 
a  specialized  area  in  which  the  literature  is  also  rather  specialized  because  of  the  complexity  of  the  fluid 
flow  problems  Though  this  subject  is  outside  the  scope  of  this  report  it  is  necessary  to  make  a  few 
brief  remarks  and  in  particular,  to  identify  the  conditions  under  which  one  must  account  for  free  sur¬ 
face  effects. 

In  the  case  of  floating  bodies  the  reader  is  referred  to  excellent  reviews  of  the  analytical  techniques 
by  W  ehausen  (Ref.  12),  Newman  (Ref. 13)  and  Oglvie  (Ref.14).  Submerged  bodies  are  only  slightly 
easier  to  handle  Some  data  on  submerged  bodies  is  given  in  Table  111.  It  should  be  stressed  that  these 
examples  are  only  pertinent  to  the  inertial  forces  generated  when  accelerating  the  bodies  ftmnwsL  Any 
prior  translation  or  rotational  motion  of  those  bodies  would  have  generated  free  surface  waves  which 
would  in  turn  affect  the  unsteady  loading  on  the  body  This  represents  the  meyor  complication  intro¬ 
duced  by  the  presence  of  a  free  surface  It  is  however  dear  that  if  the  body  motion  is  sufficiently  slow 
(characterized  by  a  velocity,  U,  say)  then  the  waves  created  will  be  negligibly  small  and  these  prior 

1 

history  effects  would  also  be  small.  This  requires  that  the  Froude  number,  UA,gd) 8  «  1. 


The  results  of  Table  Ill  do  allow  one  to  estimate  what  constitutes  proximity  to  a  free  surface  provid¬ 
ing  the  above  conditions  hold  It  can  be  seen  that  the  free  surface  has  little  effect  (less  than  5%)  pro¬ 
vided  the  ratio  of  the  depth  of  the  body  to  the  body  dimension  is  greater  than  about  4.  For  lesser 
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depths  the  added  mass  first  increases  as  the  acceleration  of  the  fluid  between  the  free  surface  and  the 
body  increases  but  then  decreases  when  the  depth  is  less  than  about  one  body  dimension  because  lesB 
fluid  is  being  accelerated 

a8  THE  EFFECT  OF  FLUID  COMPRESSIBILITY 

Generally  the  effects  of  the  conpresability  of  the  water  on  the  added  mass  can  be  neglected  in  mm* 
ocean  engineering  applications.  This  is  because  the  compressibility  does  not  begin  to  the  fluid 
flow  until  the  M  ach  number  ratio  of  the  typical  fluid  velocity  to  the  velocity  of  sound,  c,  in  the  fluid 
exceeds  a  value  of  at  least  0. 1.  In  unsteady  flows  one  must  also  consider  a  parameter  computed  as  the 
typical  acceleration  times  the  typical  body  dimenaon  and  divided  by  c.  Again  one  would  not  uorarily 
expect  any  compressibility  effect  if  this  is  less  than  0.1. 


Such  conditions  are  almost  always  satisfied  in  ocean  engineering  applications.  However  it  is  possible 
that  the  presence  of  a  large  quantity  of  bubbles  in  the  water  could  sufficiently  reduce  the  sonic  velocity, 
c,  to  such  an  extent  that  the  added  mass  would  be  altered  by  the  oompresibility  of  the  water  /gas  nix- 
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4  REVIEW  OF  EXISTING  DATA  ON  ADDED  MASS 
41  THEORETICAL  POTENTIAL  FLOW  ADDED  HASSES 

By  far  the  largest  category  of  analytical  results  for  added  mate  ere  those  calculated  for  bodies  in  tn 
infinite  fluid  domain  assuming  the  flow  to  be  potential.  The  majority  of  these  results  are  obtained  by 
methods  analogous  to  those  described  in  Section  3.  Bodies  for  which  the  steady  flows  can  be  generated 
by  superposition  of  an  array  of  potential  flow  singularities  (sources,  sinks,  doublets,  potential  vortices, 
etc)  are  particularly  compatible  with  the  use  of  expression  ( 16).  Such  methods  are  described  in  Ref-9 
and  in  many  mechanics  texts  (eg.  Ref.  1,  p.  104) .  A  particularly  useful  tabulation  of  many  of  the  avail¬ 
able  results  is  given  in  a  paper  by  Patton  (Ref. 8)  and  his  Tables  1  and  2  are  reproduced  here  as  Tables  1 
and  HI.  Note  that  the  '.ml  column  of  these  tables  contains  the  added  maas  denoted  by  the 
values  given  correspond  to  the  diagonal  terms  in  the  added  mass  matrix.  My,  the  direction  of 
acceleration  ifcnog  spec-died  in  the  second  column  (No  off-diagonal  components  of  the  added  maos 
matrix  are  Listen. .  muk  results  are  also  listed  for  bodies  on  or  near  to  a  solid  or  free  surface  and  com¬ 
ment  on  these  is  delayed  until  later.  Patton  has  included  both  theoretical  potential  flow  added  masses 
and  exj'erimentally  determined  added  masses  in  Tables  I  and  III.  These  are  distinguished  by  the  letters 
T  and  E  in  the  fourth  column  of  these  tables.  Another  excellent  source  of  tabulated  added  masses  is 
given  in  a  DTM  B  report  by  Kennard  (Ref. 9)  Kennard’s  tables  for  added  mass  coefficients  are  attached 
to  this  report  as  T  ables  11,  IV  and  V 

Though  not  exhaustive  Tables  1  through  IV  provide  a  substantial  reference  list  of  added  musses  It 
could  be  argued  with  some  justification  that  these  tables  are  more  than  adequate  for  most  engineering 
purposes  provided  the  body  under  analysis  is  not  in  dose  proximity  to  a  solid  or  free  surface.  The 
remainder  of  this  report  will  concentrate  on  the  limitations  of  this  analytical  knowledge  in  term  of 
boundary  effects  and  real  fluid  effects  (e.g.  viscous  effects).  However  before  proceed^  to  these  dis¬ 
cussions  two  further  points  should  be  made. 


First  Tables  1  through  IV  could  be  supplemented  by  the  potential  flow  methods  described  in  Section 
3  and  detailed  in  many  references  (e.g.  Ref .9).  Modem  potential  flow  computer  programs 
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(e.g.  Douglas- Neuman  code)  for  steady  flows  could  readily  be  adapted  for  this  purpose  as  discussed  in 
Section  3.5.  The  capability  to  do  this  might  be  important  in  circumstances  where  accurate  added 
masses  are  required  for  bodies  of  unusual  or  complex  geometry  or  in  circumstances  where  the  off- 
ch agonal  terms  in  the  added  mass  matrix  are  deemed  important  (the  tables  oontain  virtually  no  informer 
tion  on  off-diagonal  terms). 

The  second  point  is  that  approximate  values  for  the  conventional  or  diagonal  added  mass  terms  for 
bodies  of  complex  geometry  (for  example,  an  airplane)  can  be  obtained  by  combining  the  added  masses 
for  each  component  of  the  structure  (wings,  fualage,  tail,  etc).  Such  a  strategy  is  outlined  in  Section 

4.3. 

4.2  SENSITIVITY  TO  THE  GEOMETRY  OT  THE  BODY 

The  diagonal  terms  in  the  added  mass  matrix  (i.e  the  conventional  added  masses)  are  relatively 
independent  of  the  precise  geometry  of  a  body.  For  example,  when  accelerated  normal  to  their  longitu¬ 
dinal  axes,  cylinders  with  any  elliptical  cross-section  have  an  added  mass  equal  to  that  of  a  circular 
cylinder  with  the  same  width  normal  to  the  direction  of  acceleration  under  oonaderahon  (see  Table  I). 
Cylinders  with  more  irregular  rectangular  or  diamond  shaped  cross-sections  deviate  somewhat  from  this 
rule;  however  the  deviations  are  rather  unpredictable.  Compare  for  example  the  fact  that  the  rectangu¬ 
lar  and  diamond  shapes  in  Table  1  show  opposite  trends  as  the  cross-section  becomes  more  streamlined 
in  the  direction  of  acceleration  W  hen  the  ratio  of  cross-sectional  dimension  in  the  direction  of 
acceleration  to  that  normal  to  the  direction  of  acceleration  is  about  5  the  rectangular  shape  has 
increased  its  added  mass  by  a  factor  of  2  whereas  the  diamond  shape  has  decreased  its  added  mass  by 
a  factor  of  40%  The  unsubstantiated  opinion  of  the  author  is  that  the  experimental  values  would  show 
less  deviation  due  to  the  effects  of  flow  separation 

D  espite  these  deviations,  a  reasonable  first  approximation  to  the  translational  added  mass,  A#«,  for 
two  dimenaonal  bodies  (large  aspect  ratio  of  length  l,  to  cross-sectional  dimension,  )  would  be 
the  mass  of  a  cylinder  of  fluid  whose  diameter  is  the  same  as  the  width,  2a<,  of  the  projected  area  in 
the  direction  of  acceleration  av 
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(26) 


Consider  therefore  the  following  empirical  approximation  for  arbitrary  three  dimenticnai  bodies; 
that  the  added  masB  for  a  particular  direction  of  acceleration,  Z(,  is  given  by  the  volume  obtained  by 
rotating  the  projected  area  of  the  body  in  that  direction  about  an  axis  defined  by  the  smaller  erf  the  two 
principal  dimensions,  2a* ,  2b* ,  of  the  projected  area  For  an  elliptical  projected  area  this  would  yield 

bi>a»  (36) 


where  there  is  no  implied  summation  over  the  index  i  This  would  yield  a  reasonably  conservative 
approximation  for  the  preceding  case  of  the  cylinders.  However  it  would  substantially  overestimate  the 
aided  mass  for  a  body  like  a  sphere  which  has  a  small  aspect  ratio.  Then  the  above  estimate  would  be 
twice  the  potential  flow  value.  Perhaps  a  better  empirical  approximation  would  be 
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which  would  then  predict  both  the  cylinders  and  uie  sphere  correctly.  Testing  this  against  the  data  for  a 
prolate  ellipsoid  accelerating  "broadside  on"  (see  Table  IV)  we  find  a  value  of  the  added  mass  using 
Eq.(27)  which  is  within  7%  of  the  exact  value.  Further  improve.  ~^nts  could  dearly  be  made  but  are 
probably  of  minor  value  considering  the  other  uncertainties  discussed  below. 

43  BODIES  OF  COMPLEX  GEOMETRY 


The  result  (27)  of  the  previous  section  suggests  an  extension  for  the  purposes  of  evaluating  the  added 
mass  for  a  body  of  complex  geometry  (an  airplane) .  Though  it  would  require  further  detailed  analysis 
and  testing  it  would  not  be  unreasonable  to  suggest  that  a  complex  body  be  conadered  disassembled 
into  its  principal  component  parts  (wings,  fuselage,  tail)  and  that  the  added  masses  for  each  component 
be  evaluated  for  three  perpendicular  directions  of  acceleration  using  the  technique  outlined  in  the  previ¬ 
ous  section  Then  we  must  ask  whether  it  is  approximately  correct  to  amply  add  the  added  masses  for 
the  components  in  each  of  the  three  directions.  From  an  engineering  point  of  view  it  seems  reasonable 
to  do  this  However  it  is  very  difficult  to  give  any  quantitative  measure  of  the  error  in  such  an  estimate 
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due  to  the  interaction  of  the  components.  The  case  of  tiro  parallel  cylinders  touching  each  other  which 
is  detailed  in  Table  II,  provides  a  particular  harsh  test.  For  directions  parallel  to  and  normal  to  the 
plane  of  the  axes  of  the  cylinders  the  errors  would  be  35.5%  and  129%  respectively.  But  this  aimply 
demonstrates  that  the  two  cylinders  together  should  be  treated  as  a  single  component;  then  the  errors 
are  significantly  smaller  namely  29%  and  14%  respectively. 

M  uch  more  reasonable  tests  are  provided  by  the  winged  objects  in  Table  HI  and  we  shall,  in  particu¬ 
lar  examine  the  values  given  in  Item  3,  T  able  III.  Taking  the  individual  components  (two  flat  plates 
and  an  ellipsoid)  and  using  the  tabulated  added  masses  of  these  individual  components  in  the 
N= 0.5  one  arrives  at  a  value  of  K  of  1.293.  This  is  within  5%  of  the  actual  value  cdcuLated  namely 

1.24. 

Further  tests  would  be  needed  to  establish  confidence  limits  on  this  superposition  mathnri  but  it 
does  not  seem  unlikely  that  one  could  confidently  predict  potential  flow  added  masses  for  rampimr 
bodies  to  within  ±  30%  usmg  the  methods  outlined  above  and  empirical  formulae  such  es 
by  Eq.(27). 

<4  THE  EFFECTS  OF  A  NEARLY  SOLID  BOUNDARY 

The  discussion  in  Sections  4.2  and  4.3  was  confined  to  bodies  remote  from  a  solid  boundary.  It  is  dear 
from  the  various  examples  given  in  the  tables  that  the  presence  of  a  solid  boundary  can  r«n»  substan¬ 
tial  increase  in  the  added  mass.  This  is  due  to  the  necessary  increase  in  the  fluid  accelerations  primarily 
in  the  region  between  the  Quad  and  the  boundary.  For  example  from  Table  11  it  is  seen  that  the  aided 
mass  for  a  circular  cylinder  (radius,  o)  is  increased  by  a  factor  ae/Z h?  for  a  wail  at.  a  distance  h 
from  the  center  of  the  cylinder.  The  result  presented  is  only  approximate  and  requires  a/h  «  l.  If 
the  body  is  brought  closer  to  the  boundary  the  added  mass  increases  at  an  even  greater  rate  because 
one  develops  a  "film"  or  narrow  gap  between  the  body  and  the  wall  in  which  the  fluid  xnrelnywflnn  can 
be  very  large  indeed. 

From  a  practiced  engineering  point  of  view  there  is  a  paucity  of  data  for  these  extreme  conditions  of 
doee  proximity  of  a  body  to  a  solid  boundary.  A  simple  example  will  illustrate  some  of  tbs  dramatic 
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effects  of  tt»  proximity  of  a  solid  surface  on  the  inertial  forces  required  to  move  a  body  away  from  that 
gurfac^  Conader  the  two  dimensional  problem  of  a  flat  plate  of  width.  2a,  lying  on  an  ocean  floor. 
A  vertically  upward  force,  F,  per  unit  length  of  the  plate  is  applied  at  the  center  of  the  plate  to  lift  it 


ay  from  the  floor.  D ue  to  this  force  the  plate  has  risen  to  a  uniform  height,  h(t),  above  the  floor 
at  time  t  The  velocity  and  acceleration  of  the  plate  in  the  upward  direction  are  therefore  Ji/dt  end 
cfh/df  (see  Fig.  2). 


Figure  2. 
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This  problem  could  be  visualized  as  characteristic  of  any  fairly  flat  object  lying  on  the  ocean  floor. 
Typically  only  portions  of  the  undersurface  would  be  m  contact  with  the  ocean  floor.  However  one 
could  visualize  that  prior  to  application  of  the  force  there  is  a  typical  average  separation  distance,  ho. 
between  the  undersurface  and  the  object  Silting  up  of  the  object  could,  of  oourse,  make  h 0  very 
small.  In  any  case  some  hg  would  be  pertinent  to  the  moment  t= 0,  when  the  lift  force  is  applied. 
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W  e  concentrate  here  on  the  dynamics  of  the  body  while  the  separation,  h  is  very  small  compered 
with  the  Lateral  dimension.  2a,  of  the  object  because  it  will  be  seen  that  these  are  the  most  critical 
conditions.  Then  upward  velocity  of  the  plate,  dh.AH,  will  generate  much  Larger  horizontal  velocities, 
u  (Fig.2),  in  the  gap  than  vertical  velocities  and  hence  continuity  of  mass  in  the  gap  requires 

h  u--^-x  (28) 


and  the  momentum  equation  for  the  fluid  in  the  gap  in  the  absence  of  frictional  or  viscous  forces  yields 


1|e  +  2u|h 

p  dx  dx 


(29) 


where  p{x,t)  pressure  at  any  point,  x,  in  the  gap  Substitution  and  integration  yield  the  following 
form  for  the  pressure  distribution  in  the  gap 
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where  Pe  is  the  pressure  at  the  edges,  x-±a  W  hile  the  gap  is  small  pg  will  be  approximately  equal 
to  the  ambient  pressure,  pA ,  and  the  pressure  on  the  top-side  will  deviate  much  less  from  pA  than  the 
pressure  p  in  the  gap.  Consequently  by  integration  using  the  relation  for  p,  one  can  obtain  the 
downward  inertial  force  or  added  mass  force,  F.  imposed  by  the  fluid  on  the  plate  (per  unit  plate 
length) 
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Compare  this  with  the  known  inertial  force  on  the  plate  in  the  absence  of  the  solid  boundary  namely 
prra2  d2h/dF  where  d2hA)&  is  agam  used  to  represent  the  vertically  upward  acceleration  It  is  dear 
that  the  magnitude  of  the  former  given  by  Eq.(31)  will  typically  be  very  much  larger  than  the  latter  by 
the  large  factor,  aA  Of  course  the  preceding  analysis  ceases  to  be  valid  when  h  approaches  a  But 
it  is  dear  that  as  the  plate  is  raised  the  added  mass  per  unit  Length  begins  at  a  very  large  value  of  the 
order  of  pa3/h0  and  will  rapidly  decrease  with  h  so  that  it  asymptotically  approaches  a  value  of  the 
order  of  paz  when  h  is  of  the  order  a 
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It  is  of  interest  to  examine  briefly  the  consequences  of  this  beharior  of  the  fluid  Inertial  forces.  It, 
for  simplicity  one  neglects  the  mass  of  the  plate  itself,  then  the  upward  force  applied  to  the  plate  by  the 
cable  will  be  equal  and  opposite  to  the  fluid  inertial  force.  For  illustrative  purpose  suppose  a  constant 
upward  cable  force  is  applied  Then  integration  of  the  equation  of  motion  represented  by  Eq.(31)  if  F 
is  now  visualized  as  the  cable  force  yields  a  time  history  h.(t),  given  by 

(32) 

where  the  initial  conditions  h=h0  and  dh/dt=0  at  f=0  have  been  used.  It  is  readily  seen  that  this 
leads  to  a  kind  of  "catastrophid'  release  from  the  bottom  in  which  the  i^pwand  acceleration  increases 
with  time.  It  is  unlikely  therefore  that  a  constant  uplift  force  could  be  maintained  under  these  dr- 
cumstances.  Consequently  the  actual  initial  motion  would  be  dependent  on  other  factors  such  as  the 
cable  elasticity. 

The  author  has.  as  yet.  encountered  Utile  in  the  way  of  analysis  of  such  problems  and  suggests  this 
as  an  area  deserving  further  study  both  experimentally  and  analytically. 

4.5  VISCOUS  EFFECTS  ON  ADDED  HASS  AND  DRAG 
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The  previous  sections  of  this  chapter  have  deliberately  avoided  reference  to  a  further  complication 
caused  by  the  viscous  effects  in  the  flow  around  the  body  These  viscous  effects  on  both  the  fluid  iner¬ 
tial  and  fluid  drag  forces  have  been  the  subject  of  a  number  of  detailed  studies  as  represented  for  exam¬ 
ple  by  Refs.  10,  11  and  15.  The  essence  of  the  complication  is  th  ’  oertau.  ;*  £ifnes  of  flow  the 
viscous  processes  of  flow  separation  and  vortex  shedding  cause  radiaal  modifications  to  the  forces 
expected  on  the  baas  of  ample  addition  of  fluid  inertial  and  fluid  drag  forces.  The  latter  approximation 

is  embodied  in  what  is  known  as  M  orison's  equation  (Ref  16)  in  which  the  total  force  on  the  body.  P, 

* 

is  expressed  as 


dUj 
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where  Cy  is  a  lift  and  drag  coefficient  matrix,  and  A  is  a  typical  area  for  the  body.  This  equation  Is 
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1 


normally  quoted  for  only  one  direction  at  a  time  and  is  written  as 

—f--  JPACD[U\2  (34) 

where  Cd.  U  and  dU/dt  refer  to  a  drag  coefficient  the  velocity  and  acceleration  in  line  with  the 
force.  It  might  be  expected  that  both  M  and  CD  would  be  independent  of  the  specific  motion  under 
consideration.  However  Keulegan  and  Carpenter  (Ref.  10)  have  observed  experimentally  that,  this  was 
not  the  case  and  that  substantial  changes  in  M  and  Cp  occurred  as  the  rate  of  acceleration 
represented  by  the  period,  T,  of  their  sinusoidal  motion  was  increased  such  that  Uy  T/D  became  of 
order  one.  Here  Uy  is  the  typical  velocity  (the  peak  velocity  of  the  sinusoidal  motion)  and  D  is  the 
cylinder  or  plate  width  It  is  significant  that  all  of  their  data  was  obtained  within  a  range  of  Reynolds 
numbers,  Uy D/v  ( v  is  the  kinematic  viscosity)  between  5000  and  30,000.  Even  the  steady  flows  past 
bodies  in  this  Reynolds  number  regime  experience  substantial  unsteadiness  due  to  flow  separation  and 
vortex  shedding. 

Keulegan  and  Carpenter  found  that  the  "effect  i"e"  value  of  the  added  mass  for  cylinders  was  dose  to 
the  potential  flow  value  (pnD2/ 4  per  unit  length)  for  Uy  T/D  below  about  5  but  decreased  rapidly 
with  increasing  Uu  T/D  becoming  negative  for  a  range  of  Uu  T/D  between  10  and  20!  (Note  that 
we  have  subtracted  the  displaced  fluid  mass  from  their  results  to  get  the  true  added  mass  in  line  with 
the  discussion  of  Section  3.8 )  W  ith  further  increase  in  Uy  T/D  poative  values  similar  to  those  for 
UH  T/D  <  5  are  recovered  The  drag  coefficient,  CD,  shows  a  large  increase  for  the  same  range  of 
Uy  T/D  betweeen  10  and  20  Flat  plates  exhibited  a  different  pathological  behavior  of  the  added  mass 
and  dreg  coefficient  No  systematic  variations  with  Reynolds  number,  UyD/v.  could  be  detected 

Stop.  Ramberg  and  Ferer  (Ref.  15)  have  also  carried  out  experiments  on  sinusoidally  oscillating 
flows  except  that  the  body  rather  than  the  fluid  is  accelerated.  Their  results  do  not  agree  with  those  of 
Keulegan  and  Carpenter  For  values  of  Uy  T/D  between  1  and  12  they  found  that  the  fluid  inertial 
force  agreed  very  well  with  the  potential  flow  value  M  oreover  the  variations  in  the  effective  drag 
coefficient  could  be  accurately  predicted  by  considering  the  instantaneous  Reynolds  number  at  each 
point  during  the  cycle,  uang  some  appropriate  form  for  the  corresponding  instantaneous  drag  and 
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thereby  synthesizing  the  overall  drag  coefficient 

The  results  of  Stop,  Ramberg  end  Ferer  cannot  be  readily  reconciled  with  those  of  Keulogoa  and 
Carpenter  The  Reynolds  numbers  for  the  Stop,  Ramberg  end  Ferer  experiments  are  in  the  range 
between  230  and  40,000  and  are  therefore  similar  to  those  of  Keulegan  end  Carpenter.  It  is  quite  deer 
that  further  detailed  measurements  using  more  sophisticated  measurement  end  data  analysis  techniques 
are  needed  to  resolve  this  question  Though  it  has  little  value,  I  have  formed  the  very  tentative  opin¬ 
ion  that  the  experiments  and  data  reduction  techniques  used  by  Stop,  Ramberg  and  Ferer  are  superior 
to  those  of  Keulegan  and  Carpenter  and  therefore  I  would  place  more  confidence  in  their  results.  On 
the  other  hand  the  data  of  Keulegan  and  Carpenter  is  much  more  widely  known  and  used;  this  I  believe 
rray  be  unfortunate. 

For  the  present  it  is  necessary  for  engineering  purposes  to  be  aware  that  pathological  behavior  of  the 
fluid  inertial  forces  might  occur  for  body  motions  whose  typical  amplitude  is  greater  than  about  half  of 
the  body  dimension 

Before  leaving  this  subject  it  is  of  value  to  Second  a  few  of  the  results  of  the  experiments  ranted  out 
by  Sarpkaya  (Ref.ll)  He  oscillated  a  cylinder  in  a  direction  normal  to  the  direction  of  an  oncoming 
stream  of  fluid  and  observed  pathological  behavior  for  VT/D  (  where  V  is  now  the  steady  stream 
velocity)  between  about  3  and  10.  Furthermore,  in  one  of  the  few  experimental  measurements  of  ofT- 
dugonal  terms  in  the  force  matrix  he  observed  the  oscillations  in  the  force  on  the  body  perpendicular 
to  the  direction  of  oscillatory  motion  to  be  less  than  7%  of  the  steady  drag  in  that  direction 


-25- 


&  SUV  K  ART 

The  analytical  background  of  the  added  mass  matrix  describing  fluid  forces  due  to  aooalentton  of  the 
body  or  the  fluid  has  been  reviewed.  It  is  diown  that  the  use  of  this  ooncept  is  rigorously  justified  only 
in  the  case  of  linearly  superposable  fluid  motions  with  rigid  boundaries  In  the  context  of  ocean 
engineering  problems  this  restricts  the  analysis  to  that  of  potential  flow  and.  indeed,  almost  all  of  the 
theoretical  predictions  are  computed  from  potential  flow  analysis.  For  empirical  engineering  purposes 
the  concept  has  also  been  used  for  real  flows  with  boundary  layers,  separation  and  vortex  shedding. 

The  majority  of  potential  flow  calculations  of  added  mass  are  for  bodies  accelerated  in  an  infinite 
domain  of  incompressible,  inviscid  fluid  M  any  of  these  are  included  in  Tables  1  to  IV.  These  tables 
provide  a  substantial  reference  list  which  may  be  more  than  adequate  for  many  engineering  purposes 
provided  the  body  is  not  in  dose  proximity  to  a  solid  or  free  surface  boundary.  Furthermore  since  the 
added  mass  is  generally  rather  insensitive  to  the  detailed  geometry  of  the  body  we  have  some  prelim¬ 
inary  suggestions  as  to  how  the  added  mass  for  a  body  of  complex  geometry  might  be  estimated  As 
detailed  in  Sections  4.2  and  4.3  the  first  step  is  to  decompose  the  body  into  mayor  components.  The 
added  mass  of  each  may  then  be  estimated  for  each  direction  of  acceleration  from  the  principal  dimen¬ 
sions  (20, .  2^)  of  the  projected  area  in  that  direction  and  the  approximate  formula  (27).  The  added 
mass  for  each  component  in  each  direction  would  then  be  arithmetically  summed  1  believe  it  might  be 
posable  to  make  predictions  within  ±  30%  by  this  means.  If  better  accuracy  is  required  then  we  have 
indicated  how  modem  potential  flow  codes  (.e  g  Douglas-Neuman  code)  designed  to  calculate  steady 
flows  might  be  utilized  to  obtain  bettor  results 

D  ala  is  scarcer  for  the  cases  when  a  solid  boundary  or  free  surface  is  dose  to  the  body.  In  general 
modifications  are  required  if  such  a  boundary  is  within  four  mqjor  body  dimensons.  Free  surface 
effects  ere  quite  complex  and  are  not  covered  by  this  report  However  it  is  shown  that  a  solid  surface 
(ocean  bed)  can  have  a  dramatic  effect  particularly  when  a  body  is  being  raised  from  the  ocean  floes'. 

Data  is  also  scarce  for  the  off-diagonal  dements  of  the  added  mass  matrix.  Virtually  no  results  for 
these  elements  appear  in  the  tables  Consequently  it  is  almost  impossible  to  assess  whether  these 
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interaction  terms  ere  Important  in  practical  problems.  Again.  however,  use  of  the  aforementioned  com¬ 
puter  codes  would  permit  better  evaluation  of  the  need  to  oonader  the  off-diagonal  terms. 

The  relationship  between  the  forces  when  the  fluid  is  accelerating  past  the  body  as  opposed  to  the 
reverse  is  discussed  in  Section  3.6.  It  is  shown  that  a  relation  can  only  be  firmly  established  if  either  (i) 
superposability  is  possible  (e.g.  potential  flow)  or  (ii)  if  the  entire  previous  history  of  the  relative  velo¬ 
city  is  Identical  in  the  two  cases  Then  the  appropriate  fluid  mass  in  the  case  of  fluid  acceleration  is 
equal  to  the  added  mass  plus  the  displaced  fluid  mass. 

Finally  it  is  dear  that  viscous  effects  in  the  form  of  boundary  layer  separation  and  particularly  vortex 
shedding  could  possibly  cause  radical  departures  from  the  theoretical,  potential  flow  predictions.  The 
data  on  this  is  limited  and  contradictory.  For  the  present  one  can  only  point  out  that  pathological 
behavior  night  occur  in  certain  ranges  of  frequency  (or  typical  time  of  acceleration)  and  Reynolds 


number. 
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&  RECOH  H  ENDATIONS 

It  seems  appropriate  to  suggest  several  eras  of  engineering  importance  in  which  further  analytical. 

empirical  and  experimental  studies  would  provide  valuable  information 

A .  There  is  a  relative  paucity  of  good  experimental  data  in  the  open  literature  which  cen  be  used  to 
evaluate  the  real  fluid  effects  of  viscosity.  The  data  which  does  exist  is  often  contradictory.  Such 
experiments  are  not  easy  and  are  fraught  with  pitfalls.  However  both  measurement  techniques  and 
date  processing  methods  have  substantially  improved  in  the  last  five  years.  It  therefore  seems 
appropriate  to  suggest  further  experimental  programs  which  might  help  to  provide  some  solid  infer 
motion  that  the  engineer  could  use.  At  the  present  time  there  is  little  concrete  knowledge  which 
the  engineer  could  use  with  confidence. 

B.  The  theoretical  predictions  of  added  mass  from  potential  flow  provide  a  good  data  base  for  use  in 
estimating  the  diagonal  terms  in  the  added  mass  matrix.  This  data  base  would  be  utilized  to  pro¬ 
duce  empirical  methods  for  use  with  bodies  of  complex  geometry.  This  could  result  in  a  simple 
and  useful  computer  code  for  this  purpose. 

C.  There  are  however  very  few  calculated  values  for  the  ofl-diagonal  terms  in  added  mass  matrlcea  I 
therefore  recommend  that  in  order  to  build  up  some  data  base  for  off-diagonal  terms  and  in  order 
to  allow  more  accurate  evaluation  of  the  diagonal  terms  for  bodies  of  complex  geometry  some  of 
the  modem  potential  flow  computer  codes  (eg.  D ouglas-N euman)  be  adapted  to  evaluate  the 
entire  added  mass  matrix. 

D .  The  dramatic  effects  which  can  occur  during  separation  of  a  body  from  the  ocean  flow  drould  be 


further  investigated  both  analytically  and  experimentally. 
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TABLE  I 

(From  Reference  5) 

ADDED  (HYDRODYNAMIC)  MASSES  FOR  TWO-DIMENSIONAL 
POTENTIAL  FLOWS  ;  Reference  numbers  are  given 
under  SOURCE. 


TABLE  I  (continued) 
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TABLE  II 
(From  Reference  9) 

ADDED  MASSES  FOR  TWO-DIMENSIONAL  POTENTIAL  FLOWS 
(See  Reference  9  or  TABLE  V  for  notation) 


Circular  cylinder  in  translation  perpendicular  to  its  axis: 


a/\  U 


T j  «*—  pn  a2  V2,  as  in  Equation  [68iJ, 


m:- 


j  -  prr  a  , 


k  ■  1. 


Elliptic  cylinder  in  translation  parallel  to  an  axis,  called  the  o-axis,  either  a>b  as  shown 
or  b  >  a:  *  "" 


Tj  -  —  pn  b2  U2,  from  Equation  (841J, 

2 


hi j'  <=  prr  ab, 


k  «  b/a. 


Plane  lamina  in  translation  perpendicular  to  its  faces: 


□ 


T  ■=  —  pn  a2  V2y 
1  2 

kkt j'  -  prr  a2. 


Elliptic  cylinder  rotating  about  its  axis: 


as  in  Equation  (86b), 


Aiis  of 
Rotation 


T.  -  —  pn  (a2  -  b2)  2  «2,  as  in  Equation  (106c), 
16 

,  1  ,  ,  ,  (a2 -b2) 2 

/j  *>~pn  o6(a2  +  J2),  k 


2 ab  (a2  ♦  b2) 


Plane  lamina  rotating  about  its  central  axis: 


r; 

.5  r 

L 


Alia  of 
Rotation 


“  ~  P”  a*  u2,  as  in  Equation  [106a 'J, 
16 


kl'<* —  pn  a*. 

1  8  r 
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TABLE  II  (continued) 


Plane  lamina  rotating  about  one  edge: 


ri 


a 


Axis  of 
Rotation 

/ 

/ 

/ 


Tj  m~  pn  a*  ca2,  as  in  Equation  [106b '] 


with  /3  -  1, 

Apparent  increase  in  moment  of  inertia  9 pn  a4/8  3 

'  — .  .  ^  .  c 

Moment  of  inertia  of  fluid  displaced  by  a  3  „  a* /2  * 

cylinder  of  radius  a  rotating  as  if  rigid 
about  a  generator 


Fluid  inside  elliptic-cylindrical  shell  rotating  about  its  axis: 

,2  x2x2 


1  (a  -62)s  , 

T,m  —  pirab - ,  as  in  Equation  [105m], 


8 


a2  +  b2 


ly  »  —  pnab  (o2  +  62) 


i2-b2\2 


,a2  +  b2  j 


Fluid  inside  semicircular  cylindrical  shell  rotating  about  axis  of  the  semicircle 


Aiis  of 
Rotation 


T.m  —  ( - J  pa*  (o2,  as  in  Equation  [102e], 

4W2  2/ 


/j  =-  pa4, 


/  8  1  \ 

A;  -  21  —  -  —  J  -  0.621. 


Fluid  inside  equilateral  triangular  prism  rotating  about  its  central  axis: 


80K3f 


P  S 4  (L>2, 


i :  - 


.  - -  P  s 

1  i6  rj 


as  in  Equation  [103k], 


2 

k  e  . 
5 


Lamina  bent  in  form  of  circular  arc,  in  translation  at  angle  0  with  chord: 


i2sin20  + 


2  J 


U 2 


as  in  Equation  [78r]. 
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TABLE  II  (continued) 

Cylinder  with  contour  consisting  of  two  similar  circular  arcs;  see  Section  89. 


6 - 180/ 
degrees 

■  2c  degrees 


Cross-sectional  area  S»  — — -  [2(1-/)  rr  +  sin  2  #1 
^  s:n2  9  1  •* 


1.  Translation  parallel  to  chord  AB.  T}  _  —  pfc  SU2, 


k-  —  / - 1  —  -1. 

3  1/2  S 


1  2  ,2 

2.  Translation  perpendicular  to  chord  AB:  Tl*=—pkSU,  k 


1  \2/2  )  S 


1. 


Cylinder  with  contour  formed  by  two  simitar  parabolic  arcs  meeting  perpendicularly;  see 
Section  91(d):  ”  ”  L 


—  Phs 
1  3 


T.  --  kM'U2\ 
1  2  1 


4  K* 


1.  Translation  parallel  to  chord  AB:  k= -  1  =  0.525. 


8Af4 

2.  Translation  perpendicular  to  chord  AB:  A- - -  -  1  »  2.049. 

rr3 

Here  K  =  1.8541,  the  complete  elliptic  integral  of  modulus  V 1/2. 

Cylinder  whose  contour  is  formed  by  four  equal  semicircles: 


iWj  -  —  (2  +  rr)  ph2\  for  translation  in  any  direction 


T.~  —  kM,  U2, 
1  2  1 


k* - K2  -  1  -  1.100. 

2  +  it 


For  K,  see  the  preceding  case.  See  Section  91(e). 
Double  circular  cylinder,  each  cylinder  of  radius  a;  see  Section  90: 


,00  „ 


2  p  n  a2. 


-  Jr,2  \  rr2 

1.  Translation  parallel  to  line  of  axes  AB:  -  pn  a  V  I —  -  lj  ,  k »—  - 1 «  0.645. 

/  n2  \  rr2 

2.  Translation  perpendicular  to  line  of  ares  AB:  Tx  •  prra2V2  (— — lj  ,  A«~ -1-2. 


290. 
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TABLE  II  (continued) 


Cylinder  of  radius  a  sliding  along  fixed  plane  wall;  see  Section  90. 

Vr-'-lf')' 


Fluid u ' 
Wall 


A/j'  =■  pn  a2, 


k  - -  1  =  2.290. 

3 


Cylinder  of  rhombic  cross-section,  in  translation  along  a  diagonal;  see  Section  91(c). 

M'  =  p  s 2  sin  6 


Here  6  is  in  radians  and  T  stands  for  the  gamma  function. 
Rectangular  cylinder  in  translation  parallel  to  a  side;  see  Section  91(b)  for  references. 


M,"m  kMy-  apparent  increase  in  mass, 

A/jp  =  pnw2/i  or  M”  for  a  plane  lamina  of  width  w. 


h/w  -  0 


0.025  0.111 

1.05  1.16 


0.298  0.676 

1.29  1.42 


1.478  3.555 

1.65  2.00 


9.007  40.03 

2.50  3.50 


Circular  cylinder  with  symmetrical  fins: 


as  in  Equation  [91g], 


2 


Cylinder  of  radius  a  in  translation  and  instantaneously  coaxial  with  enclosing  fixed  cylinder 
of  radius  5: 
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TABLE  II  (continued) 


Cylinder  of  radius  a  in  translation  in  any  direction  across  axis  of  enclosing  fixed  square 
cylinder  of  side  a,  a/s  small;  see  Section  91(1), 


s 


3 


Fluid 


3 


8 


7\  *=  —  pna4U 

1  2 


2r;2  (  1  +  6.88- 


M,'-  prr  a2,  k  =  1+6.88 - , 

*  2 

8* 


Cylinder  of  radius  a  in  translation  in  any  direction  near  a  fixed  infinite  wall,  a/h  small: 


Fluid  j 
>/»vw>rv// 
Wall 


a9  in  Equation  [95g] 


prr  a 


k 


+  *••• 


2A2 


(Only  the  force  required  to  accelerate  the  cylinder  is  considered  here.) 

Cylinder  of  radius  a  moving  symmetrically  between  fixed  infinite  walls  A  apart,  a/h  rather 

small— 


Wall 
Fluid 


as  in  Equation  [46q] 


prr  a 


A- 1 


A2 


Plane  lamina  of  width  b  moving  symmetrically  between  fixed  infinite  rigid  walls  A  apart, 
b/h  rather  small:  "~™ 


Wall 

-t 


1  Fluid 


"J  CM 


rs/;r/y ///////  ■ 


as  in  Equation  [6511 


Wall 
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TABLE  1 1 1 
(From  Reference  5) 

ADDED  MASSES  FOR  THREE-DIMENSIONAL  POTENTIAL  FLOWS 
Reference  numbers  are  given  under  SOURCE. 
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TABLE  III  (continued) 


TABLE  III  (continued) 
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TABLE  III  (continued) 
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TABLE  IV 
(From  Reference  9) 

ADDED  MASSES  FOR  THREE-DIMENSIONAL  POTENTIAL  FLOWS 
(See  Reference  9  or  TABLE  V  for  notation) 


Sphere  in  translatory  motion 


T  ~  —  p  a3l/2,  as  in  Equation 


U  'm  —  np  a3. 

8  P  * 


1 

k  m - . 

2 


Sphere  moving  perpendicularly  to  infinite  rigid  plane  boundary,  a/h  small: 


'A  Fluid 


i-A 


T--pa3(  1  +  -  — 
3  \  3  A3 


—  np  o3, 

3 


+  ...y  U2,  as  in  Equation  [130a] 
with  <*-0, 

1  /  3  a3  \ 

^a“  [l+  *4*  ••••  1  • 

2  V  8  h3 


Only  the  force  required  to  accelerate  the  sphere  is  considered  here;  see  Section  130. 
Sphere  moving  parallel  to  infinite  rigid  plane  boundary,  a/h  small: 


I" 


j  ,  /  3  a3  ,  , 

T -  —  pa3  1+  —  —  +....  1  f/z,  as  in  Equation  [130a] 
3  \  16  ^3 


M  '■ —  np  a3, 

3  ¥ 

Sphere  moving  past  center  of  fixed  spherical  shell: 


with  a  -  90  deg. 


1  /  3  «3  \ 

■ —  1 1  + - +....  1 

2  V  16  A3  j 


n  .  63  +  2  a3  _ 

Tm  —  pa3  -  V2,  as  in  Equation  [I29e], 

3  i3-a3 


M  np  a3, 
3  K 


,  1  b3  ♦  2a3 

- — 

2  b3-a3 
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TABLE  IV  (continued) 

Prolate  spheroid  (or  ovary  ellipsoid),  a>b\  see  Section  137: 

Let  e  <•  eccentricity  of  sections  through  axis  of  symmetry, 


1  -ei 


In 


1  +e 
1-e 


-  2e  , 


0o- 


1  1  +  e 

- In 


1-e 


2  2  1-e 


(1)  Translation  “end  on”: 


2  ,  ,  «o 

7"  =  —  pn  ab2U2  - , 

3  2  —  a„ 


4  ,  «0 

M=  —  pn  abz,  kmk,  -  - 

3  P  ^  2-«- 


(2)  Translation  “broadside  on”: 


2  .  0o 

r=  —  par  aft2f/2  - 

3  2-/30 


4  ,  0o 

if  =*  —  pn  abl ,  fc  =  *  - 

3  2-0o 


(3)  Rotation  about  an  axis  perpendicular  to  axis  of  symmetry: 


Axis  o( 


7  =  -  *7<u2,  /=  —  pn  ab2  (a2  +  b2), 

2  X  5 


(a2  -  ft2)  (^0-«ot 


=  A ' 


(a2  +  b2)  [2(a2  -02)  -  (a2  +  b2)  (/30  -  «0)  ] 


See  Table  A 
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TABLE  IV  (continued) 


TABLE  tK 

Coefficients  of  Inertia  for  Prolate  Spheroid 


a/b 

*1 

Translation 

“end  on” 

k2 

Translation 

“broadside  on" 

k' 

Rotation  about 

Minor  Axis 

1.00 

0.500 

0 

1.50 

0.621 

0.094 

2.00 

0.702 

0.240 

2.51 

0.763 

0.367 

2.99 

0.122 

0.803 

0.465 

3.99 

0.082 

0.860 

0.608 

4.99 

0.059 

0.895 

0.701 

6.01 

0.045 

0.918 

0.764 

6.97 

0.036 

0.933 

0.805 

3.01 

0.029 

0.945 

0.840 

9.02 

0.024 

0.954 

0.865 

9.97 

0.021 

1  1 

0.883 

oo 

0 

B 

1.000 

Oblute  spheroid  (or  planetary  ellipsoid),  a  <  b,  see  Section  138,  where  b  *»  c: 
Let  e  a  eccentricity  of  sections  through  axis  of  symmetry, 

a 0  - - (e-Vl-e2  sin-1*), 

e3 

1  -  . 

ft0  - — (  Vl-e3  sin-,e  -e(l-e2)]. 
c3 

(1)  Translation  "broadside  on"  or  parallel  to  axis: 


o  a  f 

T  m  —  prr  ab 3  V3  - 


M<m  —  pit  ab3, 
3  K 


«o 

k-k,  ■»  — — 

*  2-o. 
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TABLE  TV  (continued) 

(2)  Translation  “edge  on”  or  perpendicular  to  axis: 

!  U  2  ,  ,  00 

T  =  —  pn  ab 2  (J2 


2-£0  ’ 


M=  —  pn  ab2, 


Pc 


k-k 


2  2-/3. 


r=—  a 
2 


k  =  k‘ 


/=  —  pff  a62  (a2  +  62), 
15 


(62-a2)  («0-/30) 


(a2  +  62)  [2(62-o2)  -(o2  +  62)(a0-/30)] 


See  Table  3 


TABLE  B 

Coefficients  of  Inertia  for  Oblate  Spheroid 


b/a 

*2 

Translation 

"edge  on" 

Translation 

"broadside  on” 

k' 

Rotation  about 

Equatorial  Axis 

1.00 

0.500 

0.500 

0 

1.50 

0.384 

0.803 

0.115 

2.00 

0.310 

1.118 

0.337 

2.50 

0.260 

1.428 

0.587 

3.00 

0.223 

1.742 

0.840 

4.00 

0.174 

2.379 

1.330 

5.00 

0.140 

3.000 

1.978 

6.00 

0.121 

3.642 

2.259 

7.00 

0.105 

4.279 

2.697 

8.00 

0  092 

4.915 

3.150 

9.00 

0.084 

5.549 

3.697 

10.00 

0.075 

6.183 

4.019 

oo 

0.000 

OO 

oo 
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TABLE  IV  (continued) 

Circular  disk  in  translation  perpendicular  to  its  faces: 


r--pa3  U 2,  as  in  Equation  (138o'l; 

3 

(apparent  increase  in  mass)  2 

(spherical  mass  of  fluid  of  radius  a)  n 


Circular  disk  rotating  about  a  diameter;  see  Section  138: 


Axis  of 
Rotation 


T  =  -  pa5*,2, 
45  H 


(apparent  increase  in  moment  of  inertia)  2 

- -  - 

(moment  of  inertia  of  sphere  of  fluid  of  3 
radius  a  or  8  pn  as/15) 


Elliptic  disk  of  ellipticity  e  in  translation  perpendicular  to  its  faces,  a  >  6;  References 
(240)  and  (235): 

e  =  —  V a2  -b2 7 


r*  —  p  «2bV2, 
3  E  r 


(apparent  increase  in  mass) 

4 

( —  pit  a2b  «  ellipsoidal  mass  of  fluid  with 
®  axes  a,  a,  b) 


-k”~  —  , 

E  ’ 


77/2 


E  =  /  Vl  -e2  sin2  0  dO ,  the  complete  elliptic  integral  of  the  second  kind  to  modulus  e 

for  table,  see  Peirce  (20). 

a/b  =  1  1.25  1.5  1.75  2  2.5  3  4  6  9 

k"  =  0.637  0.705  0.756  0.795  0.826  0.869  0.898  0.932  0.964  0.981 

Ellipsoid,  any  ratio  of  the  axes  a,  b,  c;  see  Section  141: 


Let  a0  =  abc 


r  d\ 

1  (a2  +  \)3/2  ( b2  +  \ ) 


(a2  +  \)3/2  (b2  +  \)i/2  (c2  +  k)1^2 


O* 

r  d  a 

/30  -  abc  / 

l  (a2  +  A)"2  (62  +  A)*2  (c2  + 


1/2 
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TABLE  IV  (continued) 


Yo 


_ d\ _ 

(a2  +  A)'^2  (b2  +  \)W  (c2  +  X)3/2 


(1)  Translation  parallel  to  the  o-axis: 


(2)  Rotation  about  the  o-axis: 


T  =  —  pn  abc 
15  H 


(62-o2)2  (y0  — / 30) 


2(b2  ~  c2)  +  (b2  +  c2)  (p0-y0) 


4  ,  ,  (^-'T(y  o“i®o) 

/'= —  prr  abc(b2  +  c2),  k  '=»  -  ■■  -  ■  . .  - 

15  2(64-c4)  +  (62  +  c2)2  (P0  —  yQ) 


For  the  expression  of  aQ,  /30,  y0  in  terms  of  elliptic  integrals,  see  N.A.C.A.  Report  210 
by  Tuckerman  (235)  or  Volume  I  of  Durand’s  Aerodynamic  Theory  (3).  Some  values  of  k  and  of 
k’,  distinguished  by  a  subscript  to  denote  the  axis  of  the  motion,  were  given  by  Zahm  (174). 


Fluid  inside  ellipsoidal  shell  rotating  about  its  o-axis,  any  relative  magnitudes  of  a,  b,  c 
(see  last  figure): 


T  = 


2 

15 


pn  abc 


(b2 -c2)2 
b2  +  c2 


as  in  Equation  [140f], 


—  pn  abc  (b2  +  c2), 
15 
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TABLE  IV  (continued) 

Solid  of  revolution  formed  by  revolving  about  its  axis  of  symmetry  the  limason  defined  by 
r  =  b  (s  +  cos  0)/(s 2  -  1)  where  b  and  s  are  constants.  The  curve  for  a  =  1  is  a  cardioid.  A 
few  values  of  k  are  given  by  Bateman  in  Reference  (240): 

s  =  1  1. I  1.2  2  3 

k  ^  0.578  0.573  0.569  0.548  0.527  0.500. 


TABLE  V 


(From  Reference  9) 

NOTATION  FOR  TABLES  II  AND  IV  :  SEE  ALSO  REFERENCE  9 


a,  b,  c 

e 

k 

In  translation, 


In  rotation. 


Radius  of  a  circle  or  semiaxis  of  an  ellipse  or  ellipsoid, 
or  half-width  or  width  of  a  lamina 

Ellipticity 

Coefficient  of  inertia,  a  dimensionless  constant 
apparent  increase  in  mass 


k  -■ 


mass  of  displaced  fluid  ’ 
2  T  2Ti 


or 


M'U 2  M[V2 


apparent  increase  in  moment  of  inertia 
moment  of  inertia  of  displaced  fluid  * 

2  T  2T\ 
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* 

O 
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V-2 


Moment  of  inertia  of  displaced  fluid  rotating  as  a  rigid  body 
about  the  assumed  axis  of  rotation 


/(  See  under  T  x 

M  '  Mass  of  fluid  displaced  by  body 


See  under  T  ^ 


T 


Kinetic  energy  of  fluid 


Ty  fy  My  Values  of  T,  l  \  M'lot  fluid  between  two  planes  parallel  to 

the  motion  and  unit  distance  apart,  in  cases  of  two- 
dimentional  motion 
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e 

P 
1 0 


Velocity  of  translation  of  body 
An  angle  in  radians 

Density  of  the  fluid,  in  dynamical  units 

Angular  velocity  of  rotation  of  a  body,  in  radians  per 
second. 


The  fluid  is  assumed  to  surround  the  body  and  to  be  of  infinite  extent  and  at  rest  at 
nfinity,  except  where  other  conditions  are  indicated.  In  regard  to  units,  see  Sections  18,  147. 


